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ABSTRACT 



Aims. Giant radio halos are mega-parsec scale synchrotron sources detected in a fraction of massive and merging galaxy clusters. 
Radio halos provide one of the most important pieces of evidence for non- thermal components in large scale structure. Statistics of 
their properties can be used to discriminate among various models for their origin. Therefore, theoretical predictions of the occurrence 
of radio halos are important as several new radio telescopes are about to begin to survey the sky at low frequencies with unprecedented 
sensitivity. 

Methods. In this paper we carry out Monte Carlo simulations to model the formation and evolution of radio halos in a cosmological 
framework. We extend previous works on the statistical properties of radio halos in the context of the turbulent re-acceleration model. 
Results. First we compute the fraction of galaxy clusters that show radio halos and derive the luminosity function of radio halos. 
Then, we derive differential and integrated number count distributions of radio halos at low radio frequencies with the main goal 
to explore the potential of the upcoming LOFAR surveys. By restricting to the case of clusters at redshifts < 0.6, we find that the 
planned LOFAR all sky survey at 120 MHz is expected to detect about 350 giant radio halos. About half of these halos have spectral 
indices larger than 1.9 and substantially brighten at lower frequencies. If detected they will allow for a confirmation that turbulence 
accelerates the emitting particles. We expect that also commissioning surveys, such as MS 3 , have the potential to detect about 60 radio 
halos in clusters of the ROSAT Brightest Cluster Sample and its extension (eBCS). These surveys will allow us to constrain how the 
rate of formation of radio halos in these clusters depends on cluster mass. 

Key words. Radiation mechanism: non-thermal - galaxies: clusters: general - radio continuum: general - X-rays: general 



1. Introduction 

Radio halos are diffuse Mpc-scale radio sources observed at the 
center of ~ 30% of massive galaxy clusters (e.g. Feretti 2005; 
Ferrari et al. 2008, for reviews). These sources emit synchrotron 
radiation due to GeV electrons diffusing through fiG magnetic 
fields and provide the most important evidence of non-thermal 
components in the intra-cluster- medium (ICM). 

Clusters hosting radio halos always show very recent or on- 
going merger events (e.g., Buote 2001; Schuecker et al 2001; 
Govoni et al. 2004; Venturi et al. 2008). This suggests a connec- 
tion between the gravitational process of cluster formation and 
the origin of these halos. Cluster mergers are expected to be the 
most important sources of non-thermal components in galaxy 
clusters. A fraction of the energy dissipated during these merg- 
ers could be channelled into amplification of the magnetic fields 
(e.g., Dolag et al. 2002; Briiggen et al. 2005; Subramanian et 
al. 2006; Ryu et al. 2008) and into the acceleration of high en- 
ergy particles via shocks and turbulence (e.g., EnBlin et al. 1998; 
Sarazin 1999; Blasi 2001; Brunetti et al. 2001, 2004; Petrosian 
2001; Miniati et al. 2001; Fujita et al. 2003; Ryu et al. 2003; 
Hoeft & Briiggen 2007; Brunetti & Lazarian 2007; Pfrommer et 
al. 2008, Brunetti et al. 2009). 
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A promising scenario proposed to explain the origin of the 
synchrotron emitting electrons in radio halos assumes that elec- 
trons are re-accelerated due to the interaction with MHD tur- 
bulence injected in the ICM in connection with cluster merg- 
ers (turbulent re-acceleration model, e.g., Brunetti et al. 2001; 
Petrosian 2001). An alternative possibility is that the emitting 
electrons are continuously injected by pp collisions in the ICM 
(secondary models; e.g., Dennison 1980; Blasi & Colafrancesco 
1999). 

In the picture of the turbulent re-acceleration scenario, the 
formation and evolution of radio halos are tightly connected with 
the dynamics and evolution of the hosting clusters. Indeed, the 
occurrence of radio halos at any redshift depends on the rate 
of cluster-cluster mergers and on the fraction of the merger en- 
ergy channelled into MHD turbulence and re-acceleration of 
high energy particles. In the last few years this has been mod- 
eled through Monte Carlo procedures (Cassano & Brunetti 2005; 
Cassano et al. 2006a) that provide predictions that can be veri- 
fied using future instruments. In this scenario radio halos have 
a relatively short lifetime (« 1 Gyr), and the fraction of galaxy 
clusters where radio halos are generated is expected to increase 
with cluster mass (or X-ray luminosity), since the energy of the 
turbulence generated during cluster mergers is expected to scale 
with the cluster thermal energy (which roughly scales as ~ M 5/3 ; 
e.g., Cassano & Brunetti 2005). It has been shown that the pre- 
dicted occurrence of radio halos as a function of the cluster mass 
(or X-ray luminosity) is in line with results obtained from a large 
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Fig. 1. A schematic representation of the theoretical synchrotron 
spectra of radio halos with different values of v s (and Vb). 
The colored regions indicate the frequency range of VLA and 
LOFAR observations. 



2. Statistical modelling of giant radio halos in 
galaxy clusters 

Turbulence generated during cluster mergers may accelerate 
relativistic particles and produce diffuse synchrotron emission 
from Mpc regions in galaxy clusters (e.g., Brunetti et al. 2008). 
Diffuse radio emission in the form of giant radio halos should be 
generated in connection with massive mergers and fade away as 
soon as turbulence is dissipated and the emitting electrons cool 
due to radiative losses. It is likely that the generation of turbu- 
lence and the acceleration of particles persist for a few crossing 
times of the cluster-core regions. This then yields a lifetime of 
about 1 Gyr. 

Since the physics of the proposed scenario is rather uncer- 
tain, we choose to model the properties of the halos and their 
cosmic evolution using a simple statistical approach. Through 
Monte Carlo calculations we take into account the main pro- 
cesses that play a role in this scenario. These include the rate of 
cluster-cluster mergers in the Universe and their mass ratios, and 
the fraction of the energy dissipated during these mergers that is 
channelled into MHD turbulence and acceleration of high energy 
particles (Cassano & Brunetti 2005; Cassano et al. 2006a). We 
refer the reader to these papers for details, here we briefly report 
the essential steps that enter into the calculations : 



observational project, the "GMRT radio halo survey" (Venturi et 
al. 2007, 2008), and with its combination with studies of nearby 
halos based on the NVSS survey (e.g., Cassano et al. 2008). 

The steep spectrum of radio halos makes these sources ideal tar- 
gets for observations at low radio frequencies suggesting that 
present radio telescopes can only detect the tip of the iceberg of 
their population (EnBlin & Rottgering 2002; Cassano et al. 2006; 
Hoeft et al. 2008). The recent discovery of the giant and ultra- 
steep spectrum radio halo in Abell 521 at low radio frequencies 
(Brunetti et al. 2008) allows a first confirmation of this conjec- 
ture and provides a glimpse of what future low frequency ra- 
dio telescopes, such as the Low Frequency Array (LOFARu and 
the Long Wavelength Array (LWA, e.g., Ellingson et al. 2009), 
might find in the upcoming years. 

LOFAR promises an impressive gain of two orders of magni- 
tude in sensitivity and angular resolution over present instru- 
ments in the frequency range 15-240 MHz, and as such will 
open up a new observational window to the Universe. In par- 
ticular, LOFAR is expected to contribute significantly to the un- 
derstanding of the origin and evolution of the relativistic matter 
and magnetic fields in galaxy clusters. 

The main focus of the present paper is to provide a theo- 
retical framework for the interpretation of future LOFAR data 
by quantifying expectations for the properties and occurrence of 
giant radio halos in the context of the turbulent re-acceleration 
scenario. In particular, in Sect. 2 we summarize the main ingre- 
dients used in the model calculations and provide an extension 
of the results of previous papers on the occurrence of radio halos 
in clusters (Sect. 2.1) and on the expected radio halo luminosity 
functions (Sect. 2.2). In Sect. 3 we derive the expected number 
counts of radio halos at 120 MHz and explore the potential of 
LOFAR surveys. Our conclusions are given in Sect.4. 

A ACDM (H = 70 Kms-'Mpc" 1 , Q m = 0.3, Q A = 0.7) 
cosmology is adopted throughout the paper. 



1 http://www.lofar.org 



i) the formation and evolution of galaxy clusters is computed 
through the extended Press & Schechter approach (1974, 
hearafter PS; Lacey & Cole 1993) that is based on the hi- 
erarchical theory of cluster formation. The PS mass func- 
tion shows a good agreement with that derived from N-body 
simulations, at least for relatively low redshifts and masses 
~ 10 14 - 10 15 /z 1 M (e.g., Springel et al. 2005), although 
it has the tendency to underestimate the number density of 
systems with mass > 10 15 h~ l M Q (e.g., Governato et al. 
1999; Bode et al. 2001; Jenkins et al. 2001). Given a present 
day mass and temperature of the parent clusters, the clus- 
ter merger history (merger trees) is obtained making use of 
Monte Carlo simulations. We simulate the formation history 
of ~ 1000 galaxy clusters with present-day masses in the 
range 2 • 10 14 - 6 • 10 15 M Q . This allows a statistical descrip- 
tion of the cosmological evolution of galaxy clusters and of 
the merging events with cosmic time; 

ii) the generation of the turbulence in the ICM is estimated for 
every merger identified in the merger trees. The resulting 
turbulence is assumed to be generated and then dissipated 
within a time-scale of the order of the cluster-cluster cross- 
ing time in that mergefl Furthermore, it is assumed that tur- 
bulence is generated in the volume swept by the subcluster 
infalling into the main cluster and that a fraction, r],, of the 
PdV work done by this subcluster goes into the excitation of 
fast magneto-acoustic waves. The PdV work is estimated as 
« pv^nr^R v , where p is the ICM density of the main cluster 
averaged over the swept cylinder, v,- is the impact velocity 
of the two subclusters, r s is the stripping radius (see also 
Sect.2.1), and R v is the virial radius of the main cluster (see 
Cassano & Brunetti 2005 for details). 

Hi) the resulting spectrum of MHD turbulence generated by the 
chain of mergers in any synthetic cluster and its evolution 
with cosmic time is computed by taking into account the in- 
jection of waves and their damping in a collisionless plasma. 



2 The cascading time scale of large scale turbulence is expected to 
be of the same order of cluster-cluster crossing time (e.g., Cassano & 
Brunetti 2005; Brunetti & Lazarian 2007). 
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Fig. 2. Ratio v s /Vb as a function of BjB cm b calculated with ho- 
mogeneous models. Particle acceleration is calculated via TTD 
resonance with magneto acoustic waves (Brunetti & Lazarian 
2007). The shadowed region marks the range of v s /vb obtained 
assuming different energy density of magneto acoustic turbu- 
lence, 3% - 30% of the thermal energy density, and duration of 
the acceleration process, At = 2-4 XT acc , where the acceleration 
time-scale is r acc = x~ X = p 2 /{4D pp ). 



Acceleration of particles by this turbulence and their evolu- 
tion is computed in connection with the evolution of syn- 
thetic clusters by solving Fokker-Planck equations and in- 
cluding the relevant energy losses. 
iv) This procedure allows for the exploration of the statistical 
properties of radio halos. Following Cassano et al. (2006), 
we consider homogeneous models, i.e., without spatial vari- 
ation of the turbulent energy, acceleration rate and magnetic 
field in the halo volume. We assume a value of the magnetic 
field, averaged over a region of radius Rh — 500 kpc, 
which scales with the virial mass of clusters, M v as: 



< B > = B < 



U M >> 



(1) 



where b > is a parameter that enters in the model calcula- 
tions. Eq.Q]is motivated by numerical cosmological (MHD) 
simulations that found a scaling of the magnetic field with 
temperature or mass of the simulated clusters {e.g., Dolag et 
al. 20020. 



2.1. Occurrence of radio halos in galaxy clusters 

Stochastic particle acceleration by MHD turbulence is believed 
to be rather inefficient in the ICM. Consequently electrons can 
be accelerated only up to energies of m e c 2 y max < several GeV, 
since at higher energies the radiation losses are efficient and 
hence dominate. The resulting gradual steepening in the theo- 
retical synchrotron spectrum of radio halos at high frequencies 
is consistent with the observed spectral shapes (or with the very 
steep spectra) of several well studied halos {e.g., Schlickeiser et 
al 1987; Thierbach et al. 2003; Brunetti et al. 2008; Dallacasa et 
al. 2009). 



3 Dolag et al. (2002) found a scaling B cc T 2 , that would imply B oc 
M 4/i if the virial scaling M cc T 3/1 is assumed. 



This steepening will make it difficult to detect these sources 
at frequencies larger than the frequency, v s , at which the steepen- 
ing becomes severe (see FigQ]). v s is expected to be a few times 
larger than the break frequency, Vb- Vb depends on the accelera- 
tion efficiency in the ICM, x, and is defined as {e.g., Cassano et 
al. 2006) : 



v b oc< 



<B> X 2 



( < B > 2 +B 2 X 

\ cmb } 



(2) 



The Transit Time Damping (TTD) is the most important colli- 
sionless resonance between the magnetosonic waves and par- 
ticles, and is due to the interaction of the compressible com- 
ponent of the magnetic field of these waves with the particles 
{e.g., Schlickeiser & Miller 1998; Cassano & Brunetti 2005; 
Brunetti & Lazarian 2007). In this case x — 4D pp /p 2 , where 
p is the momentum of the electrons and D pp is the electron dif- 
fusion coefficient in the momentum space due to the coupling 
with turbulent waves. Cassano & Brunetti (2005) derived that in 
the case of a single merger between a cluster with mass M v and 
a subcluster of mass AM, x can be approximated as: 



X ■ 



T] t I M v + AM \ 3 / 2 r 



Tjt ( M v + 



{ {Rn/r s 



if r s < R H 
s ) 1 if r s > Rh 



(3) 



where r s is the stripping radius of the subcluster crossing the 
main cluster, i.e., the distance from the center of the subcluster 
where the static pressure equals the ram pressure (see Cassano & 
Brunetti 2005 for details); Rh is the size of the radio halo, and R v 
and T are the virial radius and temperature of the main cluster, 
respectively. 

Combined with Eq. [2] this implies that larger values of v& are 
expected in the more massive clusters, oc (M v /R v ) 3 /T oc M^ 3 
(here considering for simplicity a fixed value of B, see Cassano 
et al. 2006 for more general discussion), and in connection with 
major merger events, v\, oc (1+AM/M) 3 (also r s inEq.3 increases 
with AM/M). 

Monte Carlo simulations can now be used to follow cluster- 
mergers and to explore how different mergers contribute to 
the acceleration (efficiency) of relativistic particles in the ICM. 
Consequently this allows for a statistical modeling of vj within 
a synthetic cluster sample and the derivation of its statistical de- 
pendence on cosmic time and cluster mass. 
Surveys cannot reveal radio halos that have v s smaller than the 
observing frequency, since the spectrum of these halos should 
be very steep and their emission would fall below the survey 
detection limit (Figure 1). In order to investigate the statistical 
behavior of the population of radio halos at different frequencies 
we only take halos to be observable that have v s > v . Figure 
[2] shows the ratio v s /vt calculated for homogeneous models of 
radio halos, defining v s as the frequency where the synchrotron 
spectrum of these halos is a — 1 .9 {a — 1.9 calculated between 
Vj/2.5 and v, to mimic 600-1400 MHz spectra); as Vj/vj is only 
mildly dependent on magnetic field strength and the assumed 
fraction of turbulent energy injected, we adopt v s ~ Iv/,. A sta- 
tistical modeling of v\, provides a statistical evaluation of v s in 
the synthetic cluster sample. 

In the context of the turbulent acceleration model for giant ra- 
dio halos energetics arguments imply that halos with v 5 > 1 
GHz must be generated in connection with the most energetic 
merger-events in the Universe. Only these mergers may allow for 
the efficient acceleration necessary to have relativistic electrons 
emitting at these frequencies (Cassano & Brunetti 2005). Present 
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Fig. 3. Fraction of clusters with radio halos, with v s > v a , as a function of the cluster mass in the redshift range 0-0.1 (left panel) 
and 0.4 - 0.5 (right panel). Calculations assume v„ = 1 .4 GHz, 240 MHz, 150 MHz, 120 MHz and 74 MHz (from bottom to top). 




Fig. 4. Fraction of clusters with radio halos with v s > 120 MHz (black, upper, solid lines) as a function of the cluster mass in the 
redshift range 0-0.1 (left panel) and 0.5 - 0.6 (right panel). The fractions of clusters with radio halos v s in different frequency 
ranges are also shown : v, > 1400 MHz, 600 < v s < 1400 MHz, 240 < v s < 600 MHz, and 120 < v < 240 MHz (from top to 
botton). 



surveys carried out at v ~ 1 GHz detect radio halos only in the 
most massive and merging clusters (e.g., Buote 2001, Venturi et 
al. 2008), and their occurrence has been used to constrain the 
value of r], as 0.1 - 0.3 in the models (Cassano & Brunetti 2005). 
Similar energetics arguments can be used to claim that radio 
halos with smaller values of v, must be more common, since 
they can be generated in connection with less energetic phenom- 
ena, e.g., major mergers between less massive systems or minor 
mergers in massive systems (e.g., Eqs. [2JO, that are more com- 
mon in the Universe (e.g., Cassano et al. 2008). 

In Figf3]we plot the fraction of clusters that host radio halos 
with v s > v as a function of the cluster mass and considering 
two redshift ranges : 0-0.1 (left panel) and 0.4 - 0.5 (right 
panel); this is obtained assuming a reference set of model pa- 



rameters, namely: < B >= 1.9 /J.G, b = 1.5, r\ t = 0.2 (see also 
Cassano et al. 2006). As expected, the fraction of clusters with 
halos increases at smaller values of v„, and the amount of this in- 
crement depends on the considered mass and redshift of the par- 
ent clusters, being larger at smaller cluster masses and at higher 
redshifts. 

Fig|4]we plot the fraction of radio halos with v s > 120 MHz 
(black upper line) and the differential contribution to this frac- 
tion from radio halos with v s in four frequency ranges (see fig- 
ure caption for details). For nearby systems (Fig|4] Left Panel), 
a significant fraction of massive clusters, M v > 10 15 M , is ex- 
pected to host radio halos with v s > 120 MHz; a sizeable frac- 
tion of them with v s > 600 MHz (blue and magenta lines). On 
the other hand, the majority of radio halos in clusters with mass 



R. Cassano et al.: Unveiling radio halos in galaxy clusters in the LOFAR era 



5 



M v < 10 15 M Q would have very steep spectra if observed at GHz 
frequencies, v s < 600 MHz (red line and black dot-dashed line). 
Our calculations suggest that a similar situation is expected for 
clusters at higher redshift (FigHl Right Panel). Radio halos with 
larger values of v s become much rarer with increasing redshift, 
mainly because the unavoidable inverse Compton losses at these 
redshifts limit the maximum energy of the accelerated electrons 
in these systems. At z > 0.5, only merging clusters with mass 
M v > 2 ■ 10 15 M Q have a sizeable chance to host giant radio halos 
with v s > 1.4 GHz, and an increasing contribution to the per- 
centage of radio halos at higher redshift comes from halos with 
lower v s . 



2.2. The radio halo luminosity function 

The luminosity functions of radio halos (RHLFs), i.e., the num- 
ber of halos per comoving volume and radio power, with v s > 
1.4 GHz has been derived by Cassano et al. (2006a) as: 



dN H (z) = dNniz) j dP(lA) 
dVdP(lA) dMdVj dM 



(4) 



where dN h{z) I dM dV is the theoretical mass function of radio 
halos with v s > 1.4 GHz, that is obtained combining Monte 
Carlo calculations of the fraction of clusters with halos and the 
PS mass function of clusters (e.g., Eq. 18 in Cassano et al. 2006). 

dP( \ A) I dM can be estimated from the correlation between 
the 1.4 GHz radio power, P(\A), and the mass of the parent 
clusters, that is observed for radio halos {e.g., Govoni et al. 
2001; Cassano et al. 2006a). Cassano et al. (2006) discussed the 
P( 1 A)-M v correlation in the context of the turbulent acceleration 
model and have shown that the slope is consistent with the ob- 
served value (om = 2. 9±0.4)fora well constrained region of pa- 
rameter space (B < m>, b, and rj,; Figure 7 in Cassano et al. 2006); 
model parameters adopted in the present paper, < B >- 1.9 fiG, 
b — 1.5, r\ t = 0.2, fall in this range. In particular, the value of 
the derivative dP(lA)/dM in Eq. [4] depends on the set of para- 
menters (B < m>, b) that, in the case of the reference model we are 
using in this paper, sets au = 3.3. 

To derive the RHLF at frequency v a the contribution of all 
radio halos with v s > v a should be taken into account. We first 
obtained the RHLF for halos with v s in a given frequency inter- 
val, e.g., Av si , and then combine the different contributions from 
the considered intervals Av J; : 



dN H (z) 
dVdP(v ) 



Z( 



dN H (z) 
dMdV>^ s 



x , dP(v ) y\ 
'Av.A dM V„ 



(5) 



To derive the contribution to the RHLF from radio halos with 
v s > 1.4 GHz we should calculate dP(v )/dM for these halos. 
This can be estimated from the P( 1 .4) - M v correlation assuming 
a monochromatic radio power of these halos at v : 



/ 1400 MHz \ a 

Pia(Vo,M v ) = P L4 (1.4,M,)( ) , 



(6) 



where P\a(1A,M v ) is the monochromatic radio power at 1.4 
GHz from the P(1.4) - M v correlation, and a ~ 1.3 the typ- 
ical spectral index of these halos, P(v) oc v~ a {e.g., Ferrari et 
al. 2008). 

We now consider the case of halos with v s < 1.4 GHz. The bolo- 
metric synchrotron power of radio halos is expected to scale with 
v b and B (e.g., Cassano et al. 2006) : 
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Fig. 5. Radio halo luminosity function at v o =120 MHz (black 
lines) for clusters at redshift 0-0.1 (solid thick lines) and 0.5 - 
0.6 (dashed thin lines). The contribution from halos with 120 < 
v s < 600 MHz (red lines) and with v s > 600 MHz (blue lines) 
are also shown. 



From Eqsf2]|3]it is clear that clusters with the same mass M v (and 
magnetic field B) at redshift z could have different values of v b , 
depending on the merger event responsible of the generation of 
the radio halo. Yet, Eq. [7] implies that for a fixed cluster mass 
(and consequently for a fixed value of the magnetic field) the 
synchrotron radio power emitted at the break frequency, P(v/,), 
is constant. In addition, homogeneous models, that consider an 
average value of B and y max in the halo volume, also imply that 
P(Vb)Vb P(v s )v s (Fig. |2]i. From Eq. [7] we can then derive the 
monochromatic radio power at v of halos with a given v, : 

P Ys (Vo,M v ) = P v Av s ,My)(^) a = P 1a (\A,mA-X (8) 
\v ) \v ) 

This allows the evaluation of (dP(v )/dM)& Vs . starting from 

dP(\A)ldM, and thus to derive Eq.0 

We also note that from Eqs.|6]and|8]one has : 



P Vs (Vo,M v ) = Pl.4 



(v ,My)( 



1400 MHz/ 



(9) 



P(v b )v b oc Bv b => P(v b ) oc B . 



(7) 



i.e. radio halos with synchrotron spectra steepening at smaller 
frequencies will also have monochromatic radio powers at v 
smaller than those of radio halos with larger v s . 
As a relevant example, in Fig|5]we report the expected RHLF 
at 120 MHz (black lines) for z = - 0.1 (solid thick lines) and 
z = 0.5-0.6 (dashed thick lines), where we also show the relative 
contribution from halos with v s < 600 MHz (red lines) and v s > 
600 MHz (blue lines). 

As already discussed in Cassano et al. (2006), the shape of the 
RHLF flattens at lower radio powers due to the decrease of the 
efficiency of particle acceleration in less massive clusters. We 
note that halos with v, > 600 MHz (blue lines, Fig. [5} do not 
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Fig. 6. All sky integrated RHNCs for: (1) z=0.044-0.2, obtained 
by considering a minimum mass of clusters constrained at any 
Z by the XBACS X-ray flux limit (Ebeling et al. 1996) (dashed 
black line), and by combining the above mass-constraint with 
that implied by the NVSS sensitivity (following Cassano et 
al. 2008, see their Figure 3) (solid black line); (2) z=0.2- 
0.32, obtained by considering the X-ray luminosity -range of the 
GMRT cluster sample (Venturi et al. 2007, 2008). Black points 
are the observed RHNC of giant radio halos from NVSS se- 
lected clusters in the redshift range 0.044-0.2 and re-normalized 
to account for the NVSS and XBACS sky coverage (and XBACS 
completeness). Red points are the observed RHNC of giant ra- 
dio halos in the GMRT cluster sample (with redshift z=0.2-0.32) 
and re-normalized to account for the sky coverage of the GMRT 
cluster sample. 



contribute to lower radio powers in the RHLF. This is because 
higher-frequency halos are generated in very energetic merger 
events, and must be extremely rare in smaller systems and con- 
sequently their monochromatic radio power is larger than that of 
halos with v s < 600 MHz (red lines, Fig. [5]). Finally, we note 
that with increasing redshift the RHLFs decrease due to the evo- 
lution of the cluster mass function with z and to the evolution 
with z of the fraction of galaxy clusters with radio halos (Fig. [3] 
see also Cassano et al. 2006). The evolution of the RHLF with 
Z is stronger at higher radio powers, where the dominant con- 
tribution to the RHLF comes from halos with larger v s and the 
fraction of clusters hosting these halos decreases more rapidly 
with redshift (e.g., Fig . |4j) . 

3. Number Counts of Radio Halos and LOFAR 
surveys at 120 MHz 

It has been shown that model expectations of the occurrence of 
radio halos observed at v G = 1.4 GHz are consistent with the 
fraction of radio halos with cluster mass (Cassano et al. 2008) 
and with the number counts of nearby radio halos (Cassano et 
al. 2006). 



As already discussed, in this paper we adopt a reference model 
with parameters: < B >- 1.9/iG, b — 1.5, rj, — 0.2. In Fig. [6] we 
report number counts of giant radio halos expected with these 
parameters compared with radio halo counts from the NVSS sur- 
vey at low redshift, 0.044 < z < 0.2 (Giovannini et al. 1999) 
and from the GMRT Radio Halo Survey at intermediate redshift, 
Z = 0.2 - 0.32 (Venturi et al. 2007, 2008). The latter is a pointed 
survey down to 70//Jy/beam at 610 MHz of a sample of ~ 50 
galaxy clusters extracted from the REFLEX (Bohringer et al. 
2004) and eBCS (Ebeling et al. 1998, 2000) cluster catalogs. The 
clusters have z = 0.2 - 0.4 and Z/y[0.i-2.4fci'] > 5 • 10 44 erg/s (the 
X-ray sample is complete for z < 0.32; see Cassano et al. 2008). 
All halos in the survey have 1 .4 GHz follow up. Beside the fair 
agreement between expectations and observations (see caption), 
we note that the GMRT radio halo survey is sufficiently sensi- 
tive to catch relatively faint halos and to constrain the flattening 
of the distribution of number counts of radio halos (RHNC) at 
lower fluxes. 

Encouraged by these results, in this Section we derive expected 
RHNC at 120 MHz and explore the potential of the upcoming 
LOFAR surveys. 

Because in our simplified procedure the radio power of halos 
scales with a spectral slope a = 1.3 (Eqs. [6]-[9]l and the great 
majority of halos is at z ~0.2-0.4, in the following we neglect 
K-correctionQ 

3.1. LOFAR Surveys 

LOFAR will carry out surveys between 15 MHz and 210 
MHz with unprecedented sensitivity and spatial resolution 
(e.g., Rottgering et al. 2006). Also, the unprecedented (u,v) cov- 
erage of LOFAR at short baselines maximizes the instrument 
capability to detect extended sources with low surface bright- 
ness such as radio halos. These surveys will constrain models 
for the origin of diffuse radio emission in galaxy clusters. In 
this paper we assume an observing frequency v a - 120 MHz 
where LOFAR will carry out the deepest large area radio sur- 
veys (e.g., Rottgering et al. 2006). 

The crucial step in our analysis is the estimate of the minimum 
diffuse flux from giant radio halos (integrated over a scale of ~ 
1 Mpc) that is detectable by these surveys as a function of red- 
shift. This depends on the brightness profiles of radio halos that 
is known to smoothly decrease with distance from the cluster 
center (e.g., Govoni et al. 2001). Consequently the outermost, 
low brightness, regions of halos will be difficult to detect. 
However what is important is the capability to detect the central, 
brightest, regions of radio halos in the survey images. Following 
Brunetti et al. (2007), we consider a shape of the radial profile 
of radio halos that is obtained from the analysis of well studied 
halos. We assume a circular observing beam = 25x25 arcsec, 
and follow two complementary approaches^ 

i) since radio halos emit about half of their total radio flux 
within their half radius (Brunetti et al. 2007), we estimate 
the minimum flux of a detectable halo, / min (z), through the 

4 For simplicity we also consider observables those halos with v s ~ 
v„ regardless of their redshift. This would slightly affect only the num- 
ber counts of halos with v, < v a ( 1 + z) that represent a minimal fraction 
of our halo population. 

5 The 120 MHz LOFAR survey will have a full resolution of 5-6 arc- 
sec, thus we are considering the case of tapered images that increase 
the sensitivity to extended emission without changing significantly the 
point source sensitivity (due to the large number of inner LOFAR sta- 
tions). 
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Fig. 7. Minimum flux (left panel) and power (right panel) of detectable radio halos at 120 MHz using Eq.[l0](solid lines) and Eq.fTTI 
(dashed lines). %F = 0.1, 0.25, 0.5, 1.0, 1.5 mJy/b are assumed (from bottom to top), and beam=25x25 arcsec. 
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Fig. 8. Number (all-sky) of radio halos with v s > 120 MHz (black lines) as a function of redshift that can be expected at the 
sensitivity of LOFAR surveys. Calculations are performed following approaches i) in Sect. 2.1 (left) and ii) (right), assuming 
£ rms- 0.1, 0.25, 0.5, 1.0, 1.5 (bottom to top). Red lines give the number counts of radio halos with 120 < v s < 600 MHz. 



requirement that the mean brightness within Rh/2, J^jif 
is £ times the rms, F, of the survey, i.e. : 



fminiz) - 1(T 3 ( f4 Rfe) [mJy] 

\0.5mJy/beam/ 



5mJy /beam^ 



(10) 



where 8h(z) is the angular size of radio halos, in arcseconds, 
at a given redshift; allowing for the detection of diffuse halo 
emission in the images produced by the survey. Injection of 
fake radio halos in the (u,v) plane of interferometric data 
from NVSS observations show that radio halos at z < 0.3 be- 
come visible in the images as soon as their flux approaches 
that obtained by Eq.[10]with f ~ 1 - 2 (Cassano et al. 2008); 
ii) following a second approach we estimate the minimum flux 
of a detectable halo through the requirement that the average 



brightness within 5 observing beams is 3 x £ times the rms, 
F, of the survey. The minimum flux is obtained through the 
condition : 



r 

Jo 



2n I(b)bdb = 5S 

beam Off) 



(ID 



where 1(b) is the typical radial profile of halos (Brunetti et 
al. 2007), Sbeam is the beam area, and b s = (5 S beainlri) 1 ^ 2 ■ 
The aim of this second approach is to avoid any bias due to 
the redshift of the halos and is motivated by the fact that in 
the first approach the sensitivity limit is reached in a fairly 
large area (many beams) for nearby radio halos, but only in 
a few beam area in the case of halos at z=0.5-0.6. 



8 



R. Cassano et al.: Unveiling radio halos in galaxy clusters in the LOFAR era 




ct[120-300 MHz] 



Fig. 9. Spectral index distribution of radio halos, calculated be- 
tween 120-330 MHz, for the number counts in Fig.[8](left panel) 
and assuming £ F= 0. 1 (filled symbols) and 0.5 mJy/b (empty 
symbols). Spectral indices are calculated for homogeneous mod- 
els, with parameters given in Sect. 2, assuming the Vb distribu- 
tion derived from the Monte Carlo-synthetic clusters. Radio ha- 
los are binned according to the following values of Vb'. Vj, > 90 
MHz, 35 < v h < 90 MHz, and v b < 35 MHz (from right to 
left). The horizontal bars encompasse the variation in a due to 
the range of vjvb in Fig. [2] 



Fig|7] shows /„„■„ of radio halos as a function of redshift 
(left panel), and the corresponding minimum radio power (right 
panel), obtained following the two approaches and assuming 
f F= 0.1 . . .1 .5 mJy/beam (see Figure caption). 
Given the RHLF (dNn(z) /dP(v„)dV) the number counts of radio 
halos with / > / m ,„(z) in a redshift interval, Az = z% — Zu is given 
by: 



Fig. 10. Minimum X-ray luminosity of clusters with detectable 
halos at v o =120 MHz. Calculations are shown for ^ F-0.25 
(dashed lines) and 1 mJy/b (solid lines), and for different ranges 
of v s : 120 < v s < 240 MHz (black lines), 240 < v s < 600 MHz 
(red lines), 600 < v s < 1400 MHz (blue lines) and v s > 1.4 GHz 
(magenta lines) (from top to bottom). 
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In Fig. [8] we show the all-sky number of radio halos with 
v s > 120 MHz in different redshift intervals detectable by typ- 
ical LOFAR surveys with different sensitivities (0.1 ... 1.5 
mJy/beam, see figure caption) and following the approach i) 
(Left Panel) and ii) (Right Panel). 

The LOFAR all sky survey (e.g., Rottgering 2009, priv. comm.) 
is expected to reach an rms=0.1 mJy/beam at 120 MHz. 
Considering the case i) (Fig. [8] Left Panel) with £ ~ 2-3, 
we predict that this survey will detect more than 350 radio ha- 
los at redshift <0.6, in the northern hemisphere (6 > 0) and at 
high Galactic latitudes (\b\ > 20). This will increase the statistics 
of radio halos by about a factor of 20 with respect to that pro- 
duced with the NVSS. The LOFAR commissioning MS 3 survey 
is expected to reach sensitivities a;0.5 mJy/b at 150 MHz. Based 
on our results, »T00 radio halos are expected to be discovered 
through this survey within one year time-scale. 

The spectral properties of the population of radio halos visi- 
ble by the future radio surveys at low frequencies are expected to 
change with increasing the sensitivity of these surveys. In Fig. [8] 
we show the total number of halos with v s > 120 MHz (solid 



lines) together with the number of halos with a spectral steepen- 
ing at low frequencies, 120 < v s < 600 MHz. The latter class 
of radio halos has a synchrotron spectral index a > 1 .9 between 
250-600 MHz, and would become visible only at low frequen- 
cies, v < 600 MHz. We find that about 55% of radio halos in 
the LOFAR all sky survey at 120 MHz is expected to belong 
to this class of ultra-steep spectrum radio halos, while radio ha- 
los with larger v s are expected to dominate the population in 
shallower surveys. This is simply due to the fact that, for the 
reasons explained in Sect. 2.2, low frequency radio halos are ex- 
pected to populate the low power-end of the RHLF (e.g. Fig. [5j. 
<iP(v )(12S° m Pl ementar y information is given in Fig. [9] that shows the 
expected distribution of halo spectral indices, that refer to the 
number distributions in Fig. [8] and its evolution with sensitivity 
of radio observations; spectra in Fig. [9] are calculated between 
120-300 MHz assuming homogeneous models. 
Ultra-steep spectrum halos are a unique prediction of the turbu- 
lent re-acceleration model (e.g., Brunetti et al. 2008) and our ex- 
pectations demonstrate the potential of LOFAR in constraining 
present models for the origin of radio halos. 



3.2. Application to X-ray selected cluster samples 

Although unbiased surveys of radio halos provide an important 
way to address the occurrence of these sources (Sect. 3.1), a po- 
tential problem with these approaches is the identification of ra- 
dio halos and of their hosting clusters. This is because radio ha- 
los constitute only a very small fraction of the entire radio source 
population and need to be distinguished from confused regions 
of superpositions for radio AGNs and starbust galaxies. 

Alternatively, an efficient approach is to exploit deep LOFAR 
surveys at low radio frequencies of X-ray selected samples of 
galaxy clusters. 
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Fig. 11. Integrated (left) and differential (right) number counts of radio halos from radio follow up of eBCS and MACS clusters 
(see text). Calculations are shown for £ F- 0.25 (upper panels) and 1.0 mJy/b (lower panels) at 120 MHz. Thick (black) solid 
lines give the case v s > 120 MHz, while differential contributions are shown with different colors : v s > 1.4 GHz (magenta lines), 
600 < v s < 1400 MHz (blue lines), 240 < v s < 600 MHz (red lines) and 120 < v s < 240 MHz (black thin lines). 



Here we derive the number of radio halos, and their flux and 
redshift distribution, that are expected from LOFAR observa- 
tions of sample of X-ray selected clusters. 

There are several catalogs of X-ray selected clusters in 
the northern hemisphere that contain clusters extracted from 
the ROSAT All-Sky Survey (RASS, Trumper 1993). At red- 
shift < 0.3 the ROSAT Brightest Cluster Sample and its exten- 
sion at lower X-ray fluxes (eBCS, Ebeling et al. 1998, 2000) 
and the Northern ROSAT All-Sky (NORAS) Cluster Survey 
(Bohringer et al. 2000) provide cluster catalogs with X-ray flux 
/x[0.i-2.4]kev > 3 ■ 10~ 12 erg s^'cirr 2 ; the eBCS is 75% complete 
down to this flux limit. The extension towards larger redshifts of 
these catalogs is given by the Massive Cluster Survey (MACS, 
Ebeling et al. 2001) that contains clusters with /A-[o.i-2.4]keV > 
1 • 10~ 12 erg s"'cm" 2 at z=0.3-0.6. All these surveys have optical 



follow-ups and provide a useful starting point to pick up radio 
halos in LOFAR surveys. 

A well-known correlation exists between the synchrotron 
power of giant radio halos and the X-ray luminosity of the 
hosting clusters, P(1.4) oc L x x , with x 2 (e.g., Liang et 
al. 2000; Bacchi et al. 2003; EnBlin & Rottgering 2002; Cassano 
et al. 2006; Brunetti et al. 2009). This implies that the X-ray flux 
limit of the survey, fx, is related to the radio flux of halos. The 
minimum flux of radio halos that can be detected at redshift z is 
given by the maximum value between the minimum radio flux 
due to the sensitivity of radio surveys (Sect. 3.1) and that con- 
strained by fx through the radio - X-ray correlation. To address 
this issue at v o =120 MHz in the case of radio halos with v s > 1 .4 
GHz we assume a correlation between the monochromatic radio 
luminosity at 120 MHz and Lx rescaled from that at 1.4 GHz 
through Eq. [6] For halos with lower v s (yet v s > 120 MHz) the 
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Fig. 12. Number distributions of radio halos in eBCS clusters from radio follow up are shown in different X-ray luminosity bins 
of the hosting clusters and in two redshift bins : 0-0.2 (left panel) and 0.2-0.3 (right panel). Calculations are obtained assuming a 
sensitivity at 120 MHz, f F— 1 mJy/b and beam = 25x25 arcsec. Solid lines show the number of radio halos with v s > 120 MHz, 
while the dashed areas mark the contribution from halos with 120 < v s < 600 MHz. For comparison, the dashed lines show the 
expectations assuming that a fixed fraction of galaxy clusters = 30% host radio halos and that all radio halos have v s > 1400 MHz. 



correlation between the radio luminosity at 120 MHz and Lx is 
obtained via Eq.[8]that accounts for the smaller radio power ex- 
pected in the case of halos with steeper spectrum (Sect. 2.2). 

In this Section we model the sensitivity of LOFAR at 120 MHz 
following the approach i) in Sect. 3 (Eq.fTOb. More specifically, 
to detect radio halos, we consider 120 MHz LOFAR follow ups 
of a cluster catalog obtained by combining the eBCS (at z < 0.3) 
and the MACS (0.3 < z < 0.6) samples, and assume references 
sensitivities of radio observations £ F = 0.25 and 1 mJy/beam. 
The minimum Lx of clusters where these radio observations are 
expected to detect giant radio halos results from the combina- 
tion of the above radio sensitivity and the minimum Lx in clus- 
ter catalogs at redshift z, and is shown in Fig.[10]by considering 
different v s (see Figure caption for details). One may note that 
at intermediate redshift and at higher redshift the luminosity- 
limit is driven by the X-ray flux limit of the eBCS and MACS 
catalogs, respectively. On the other hand, we expect that in the 
redshift range where the minimum Lx is constrained by the ra- 
dio sensitivity, radio halos with v s in the range 120-240 MHz 
(Fig. [T0l black lines) can be detected in clusters with X-ray lu- 
minosity about 3 times larger than that of clusters with v s > 1 .4 
GHz halos (Fig. [TO] magenta lines). 

In FigQ~TJwe show the cumulative and differential number counts 
of radio halos expected from the LOFAR follow up of eBCS and 
MACS clusters at 120 MHz. This is obtained through Eq. [12] 
and by taking into account both the selection criteria shown in 
Fig[10]and the sky coverage of eBCS and MACS surveys. The 
inflection in the number counts at z — 0.3 is due to the change in 
the X-ray selection criteria (see FigfTOb moving from the eBCS 
(z < 0.3) to the MACS (z > 0.3) cluster sample. We expect 
that the LOFAR all sky survey, with a planned sensitivity in line 
with the case £ F — 0.25 mJy/beam (FigQT| upper panels), will 
discover about 130 radio halos out of the ~ 400 clusters in the 
eBCS and MACS catalogs. Remarkably, about 40% of these ra- 
dio halos are expected with v, < 600 MHz, thus halos with ex- 



treme steep spectra at GHz frequencies. The majority of radio 
halos in eBCS and MACS clusters is expected at z = 0.2-0.4, 
while the small number of clusters at z > 0.5 with X-ray flux 
above the flux limit of the MACS catalog does not allow a sta- 
tistically solid expectations, although we may expect a couple 
of radio halos hosted in MACS clusters at this redshift. At this 
redshift we expect that only major mergers in massive clusters 
(M„ > 2 • 10 15 M ) can generate radio halos with v s > 1.4 GHz 
(Fig.HJ right panel). The powerful radio halo recently discovered 
in the cluster MACS J0717.5 +3745 (e.g., Bonafede et al. 2009; 
van Weeren et al. 2009) is consistent with these expectations. 

Fig. QT| (lower panels) shows the expected number counts of ra- 
dio halos assuming the more conservative case £F = 1 mJy/b 
that is suitable to explore the potential of the LOFAR MS 3 com- 
missioning survey. In this case about 80 radio halos are expected 
to be found in eBCS and MACS clusters; about 20 of these halos 
are expected with v, < 600 MHz. We notice that the number of 
expected radio halos from follow up of eBCS and MACS clus- 
ters increases by less than a factor of 2 considering a substantial 
drop in radio sensitivity from %F — 1 to 0.25 mJy/b. This is not 
surprising as the majority of radio halos that are expected to be 
discovered by deep radio observations should be found in galaxy 
clusters with X-ray luminosity below the luminosity-threshold 
of the eBCS and MACS catalogs {e.g., Fig. [TOt . 

The eBCS cluster sample contains 300 galaxy clusters at 
z < 0.3 and covers the northern hemisphere. The redshift and X- 
ray luminosity distribution of eBCS clusters is public (Ebeling 
et al. 1998, 2000) and thus we can provide a more quantitative 
expectation based on e.g., the more conservative case, MS 3 -like 
that assumes %F = 1 mJy/beam at 120 MHz (in this case the 
selection function of clusters in the Lx-z plane is reported in 
Fig. [10] solid lines at z < 0.3). 

In FigQ~2]we show the distribution of the expected radio halos in 
the eBCS clusters in two redshift intervals: - 0.2 and 0.2 - 0.3 
(left and right panels, respectively). We find that radio observa- 
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tions at 120 MHz are expected to discover radio halos in about 60 
clusters, i.e., in about 20% of eBCS clusters. In addition, about 
12 of these halos are expected to have very steep radio spectra, 
v, < 600 MHz (magenta shadowed region in Fig.[T2l. 
Finally, the percentage of clusters with radio halos is expected 
to increase with the X-ray luminosity of the hosting clusters. 
This is particularly relevant in the redshift interval z = - 0.2 
when compared to expectations calculated under the assumption 
that the fraction of clusters hosting radio halos is constant with 
cluster mass (FigQ~2] dashed lines, see caption). Consequently 
LOFAR will be able to readily test this unique expectation of the 
turbulent re-acceleration model. 

4. Summary & Conclusions 

In the present paper we perform Monte Carlo simulations to 
model the formation and evolution of giant radio halos in the 
framework of the merger-induced particle acceleration scenario 
(see Sec. 2). Following Cassano et al. (2006a) we use homo- 
geneous models that assume a) an average value of the mag- 
netic field strength in the radio halo volume that scales with 
cluster mass as B = B < M>M b v , and b) that a fraction, 77,, of the 
PdV work, done by subclusters crossing the main clusters dur- 
ing mergers goes into magneto-acoustic turbulence. Although 
simple, these models reproduce the presently observed fraction 
of galaxy clusters with radio halos and the scalings between 
the monochromatic radio power of halos at 1.4 GHz and the 
mass and X-ray luminosity of the host clusters (e.g., Cassano 
et al. 2006, 2008; Venturi et al. 2008), provided that the model 
parameters (B < m>, b, T] t ) lie within a fairly constrained range of 
values (Fig. 7 in Cassano et al. 2006a); in the present paper we 
adopt a reference set of parameters: < B >= 1.9/iG, b = 1.5, 
rjt = 0.2, that falls in that range. 

In Fig. |6]we show that the expected number counts of giant 
radio halos at v„ = 1.4 GHz obtained with this set of param- 
eters are nicely in agreement with both the data at low redshift 
(NVSS-XBACS selected radio halos, Giovannini et al. 1999) and 
at intermediate redshift (clusters in the "GMRT radio halo sur- 
vey", Venturi et al. 2007, 2008). 

The most important expectation of the turbulent re- 
acceleration scenario is that the synchrotron spectrum of radio 
halos should become gradually steeper above a frequency, v s that 
is determined by the energetics of the merger events that gener- 
ate the halos and by the electron radiative losses (e.g., Fujita et 
al. 2003; Cassano & Brunetti 2005). Consequently, the popula- 
tion of radio halos is expected to be constituted by a mixture 
of halos with different spectra, with steep spectrum halos be- 
ing more common in the Universe than those with flatter spectra 
(e.g., Cassano et al. 2006). The discovery of these very steep 
spectrum halos will allow to test the above theoretical conjec- 
tures. 

In Sect. 2 we have derived the expected radio halo luminosity 
functions (RHLF) at frequency v , that account for the contribu- 
tion from the different populations of radio halos with v s > v a . 
The RHLF are obtained combining the theoretical mass function 
of radio halos (with different v s > v a ) with the radio power- 
cluster mass correlation (EqHJi. The expected monochromatic 
radio power at v of halos hosted by clusters with mass M v is 
extrapolated from the observed P( 1 A)-M v correlation by assum- 
ing simple scaling relations, appropriate for homogenous mod- 
els, that account for the dependence of the emitted synchrotron 
power on v s (Eqs |8l9l ). 

As a relevant case we calculate the expected RHLF at v = 120 
MHz (Fig. |5). The shape of the RHLF can be approximated by 



a power law over more than two orders of magnitude in radio 
power. Homogeneous models imply the following scalings be- 
tween v s , cluster mass and the radio luminosity at v , P Vs (y<>) '■ 



M 4/3+fe. 



(1 + AM/My 



«b> 2 +Bi mh r 

and from Eq.|9]and the Pi 4 - M v correlation: 

p v ,(v ) oc My s 



(13) 



(14) 



i.e., radio halos with larger v s are typically generated in massive 
clusters that undergo major mergers and contribute to the RHLF 
at larger powers. On the other hand halos with smaller v s are 
typically generated in less massive systems and contribute to the 
RHLF at fainter powers. Radio halos with v s > 120 MHz how- 
ever become increasingly rare in clusters with mass < 5xl0 14 M G 
and this explains the drop of the RHLF at lower radio powers in 
Fig- El At the same time, halos with monochromatic radio emis- 
sion at 120 MHz > 10 26 W Hz 1 would be generated in connec- 
tion with very energetics merging events in very massive clus- 
ters, that are extremely rare, and this explains the RHLF cut-off 
at higher synchrotron powers in Fig. [5] 

In Sect. 3 we discuss the expected number counts of radio 
halos at 120 MHz, that best allow us to explore the potential of 
incoming LOFAR surveys in constraining present models. 
The crucial step in this analysis is the estimate of the minimum 
diffuse flux from giant radio halos that is detectable by these sur- 
veys. Because the LOFAR capabilities will become more clear 
during the incoming commissioning phase we exploit two com- 
plementary approaches: i) we required that at least half of the 
radio halo emission is above a fixed brightness-threshold, % F 
(F being the rms of LOFAR surveys; ii) we required that the 
signal from the radio halo is > 3 x £ F in at least 5 beam area of 
LOFAR observations. In both cases we assume a fixed shape of 
the radial profile of radio halos, calibrated through that of several 
well studied halos at 1 .4 GHz, that introduces a potential source 
of uncertainty. 

Although the uncertainties due to the unavoidable simplifi- 
cations in our calculations, the expected number counts of ra- 
dio halos highlights the potential of future LOFAR surveys. By 
assuming the expected sensitivity of the LOFAR all sky sur- 
vey (e.g., Rottgering 2009; priv. comm.), rms =0.1 mJy/b, and 
£ ~ 2 - 3, we predict that about 350 giant radio halos (~ 200 
considering the case ii)) can be detected at redshift <0.6. This 
means that LOFAR will increase the statistics of these sources 
by a factor of ~ 20 with respect to present day surveys. About 
55% of these halos are predicted with a synchrotron spectral in- 
dex a > 1.9 between 250-600 MHz, and would brighten only 
at lower frequencies, unaccessible to present observations. Most 
important, the spectral properties of the population of radio halos 
are expected to change with increasing the sensitivity of the sur- 
veys as steep spectrum radio halos are expected to populate the 
low-power end of the RHLF. A large fraction of radio halos with 
spectrum steeper than a x 1.5 (e.g., Fig. |9]l is expected to allow 
a prompt discrimination between different models for the origin 
of radio halos, for instance in this case simple energetic argu- 
ments would rule out a secondary origin of the emitting electrons 
(e.g., Brunetti 2004; Brunetti et al. 2008). 

Due to the large number of expected radio halos, a potential 
problem with these surveys is the identification of halos and of 
their hosting clusters. As a matter of fact we expect that LOFAR 
surveys will open the possibility to unveil radio halos in galaxy 
clusters with masses > 6 - 7 x 10 14 M G at intermediate red- 
shift. On the other hand, statistical samples of X-ray selected 
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clusters, that are unique tools for the identification of the hosting 
clusters, typically select more massive clusters at intermediate 
z. Consequently, we explored the potential of the first LOFAR 
surveys as deep follow ups of available X-ray selected samples 
of galaxy clusters. 

We calculate the radio halo number counts expected from the 
follow up of clusters in the eBCS and MACS samples that collect 
~ 400 galaxy clusters in the redshift range 0-0.6. We expect 
that the LOFAR all sky survey, with a planned sensitivity in line 
with the case ^ F = 0.25 mJy/b, will discover about 130 radio 
halos in eBCS and MACS clusters and that about 40% of these 
radio halos should have very steep spectrum, v s < 600 MHz. The 
majority of radio halos in eBCS and MACS clusters are expected 
at z = 0.2-0.4, while the small number of clusters at z > 0.5 in the 
MACS catalog does not allow a statistically solid expectations, 
although we may expect a couple of radio halos hosted in MACS 
clusters at this redshift. 

The MS 3 survey will be carried out in 2010, covering the 
northern hemisphere, and is expected to reach a noise level 
of about 0.5 mJy/b at 150 MHz, that implies a sensitivity 
to diffuse emission from galaxy clusters about one order of 
magnitude (assuming a as 1.3) better than present surveys 
(e.g., NVSS, Condon et al. 1998; VLSS, Cohen et al. 2007; 
WENSS, Rengelink et al. 1997). 

We considered MS 3 pointings relative to the fields of the about 
300 galaxy clusters at z < 0.3 in the eBCS catalogues. We find 
that about 60 radio halos are expected to be detected by MS 3 
observations in these clusters, 25% of them (10-15 halos) are 
expected with v s < 600 MHz. Fairly sensitive GMRT observa- 
tions of eBCS clusters at redshift 0.2-0.3 are already available 
(Venturi et al. 2007, 2008) and we expect that in a few cases ra- 
dio halos would glow up in the MS 3 images where no diffuse 
radio emission is detected at 610 MHz. We also find that MS 3 
observations of eBCS clusters at z = - 0.2 can be used to test 
the increase of the fraction of cluster with radio halos with the X- 
ray luminosity of the host clusters, which is a unique prediction 
of our model (FigfT2l 

The most important simplification in our calculations is the 
use of homogeneous models. Non-homogeneous approaches, 
that model the spatial dependence of the acceleration efficiency 
and magnetic field in the halo volume (e.g., Brunetti et al. 2004), 
and possibly their combination with future numerical simula- 
tions, will provide a further step to interpret LOFAR data. Also 
the use of the extended PS theory is expected to introduce some 
biases. For instance, it is well-known that the PS mass func- 
tion underproduces the expected number of massive clusters 
(M > 10 15 M ) at higher redshift, z ~ 0.4-0.5, by a factor of ~ 2 
with respect to that found in N-body simulations (e.g., Governato 
et al. 1999; Bode et al. 2001; Jenkins et al. 2001). Since in our 
model the great majority of halos at these redshift is associated 
with massive clusters, the use of the PS mass function implies 
that the RHNC at z > 0.4 - 0.5 could be underestimated by 
a similar factor. A further refinement of the approach proposed 
in the present paper could be obtained with the use of galaxy 
cluster merger trees extracted from N-body simulations. These 
would also allow for a more realistic description of the merger 
events (spatially resolved, multiple mergers, etc.). 

In the present paper we focus on a reference set of model 
parameters. Cassano et al. (2006a) discussed the dependence of 
model expectations at 1.4 GHz on these parameters. Based on 
their analysis we expect that all the general results given in the 
present paper are independent of the adopted values for param- 
eters. The expected number counts of halos should change by a 
factor of ~ 2 - 2.5 considering sets of model parameters within 



the region (< B >,b,i] t ) that allows for reproducing the observed 
Pi a - M v correlation. In this case the number of halos that we 
expect decreases from super-linear sets of parameters (b > 1 and 
< B >> 1.5 fiG) to sub-linear cases (b < 1 and < B >< 1.5 fiG) 
(see also Fig.4 in Cassano et al. 2006b); a more detailed study 
will be presented in a future paper. 
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